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Abstract 

Dynamics of nanoscopic arrays of monodomain magnetic elements is 
simulated by means of the Pardavi-Horvath algorithm. Experimental hys- 
teresis loop is reproduced for the arrays of Ni, with the period 100 nm 
and the mean coercive field 710 Oe.We investigate the box-counting frac- 
tal dimension of a cluster of elements with given orientation of magnetic 
moments. No fractal behavior is found. Also, the damage spreading tech- 
nique is applied to check the criticality. We find that the consequences 
of a local flip of one magnetic element remainlimited to a finite area. We 
conclude that the system does not show a critical behavior. 

1 Introduction 

Nanoscopic magnetic arrays have been proposed recently as devices of ultrahigh- 
density information storage 1[. On the other hand, it was announced recently 
OH that an array of bistable magnetic wires can display a complex behavior. 
Namely, the authors discussed the effect of self-organized criticality (commonly 
abbreviated as SOC) 0]. The problem is that a system in the critical state 
cannot be used to store any information, because a flip of one magnetic moment 
can change the magnetic structure of the whole array. The aim of this work is 
to check by means of a computer simulation if the array is indeed in the critical 
state. The values of the simulation parameters are selected as close as possible 
to the experimental system described in £Q. 

The dynamics of the total magnetization, i.e. its time dependence in the 
presence of the oscillating magnetic field, is governed by two agents: the long- 
range magnetostaticinteraction between each two elements of the array, and the 
switching field(coercive field), which varies from one element to another. The 
nonzcroswitching field enables to preserve information in the array at positive 
temperatures. The interaction leads to flips of magnetic moments and thcinfor- 
mation — otherwise frozen in the array — is destroyed, at leastpartially. 

Let us recall the characteristic features of SOC. According to an illuminating 
paperby Flyvbjerg |o], excitations wandering in a system can lead it to a self- 
organized state, i.e. a spontaneously formed state far from thermalequilibrium. 
By criticality we mean that in this state, excitations can be of any order of 
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magnitude. In fact there is no characteristic scale in a criticalstate: the array 
is scale-free in the same way as a cluster of spins "up" in a ferromagnet at its 
Curie temperature. 

Basically, we search for excitations which spread over the lattice. The tech- 
nique applied is known as damage spreading [7]. It is believed, that if SOC is 
present in a system, a local variation of its structure can lead to a global change 

Then, existence of such a spreading is an argument for the presence of crit- 
icality. On the other hand, one of the arguments of Refs. |21 E] is based on the 
calculation of the box-counting fractal dimension 9 of the cluster of elements 
with a given orientation of magnetic moments. That is why we investigate also 
this fractal dimension. 

2 The system and the Pardavi-Horvath algo- 
rithm 

Our simulations concentrate on the nanoscopic array described in Ref. pQ. Mon- 
odomain magnetic cylinders of nickel, 57 nm of diameter and 115 nm length, 
form a square array with the period 100 nm. The magnetization is M s = 370 
emu/cm 3 . The mean switching field of one element is H s = 710 Oe, with the 
standard deviation 105Oe. The switching field is constant in time for each ele- 
ment. The magnetostatic interaction is calculated by means of the rectangular 
prism approximation formula 6 . The only modification introduced to the in- 
vestigation of SOC are the periodic boundary conditions, to avoid boundaries 
which could alter the results. 

The dynamics of the system is simulated by means of the Pardavi-Horvath 
(PH) algorithm 10^. Initial parameters of this algorithm are an external mag- 
netic field H ex , the magnetic states of each element of the array and their 
switching fields H s . At the begining usually H ex — and the magnetic states 
of the elements are random. Then for every element of the array, the interac- 
tion field Hi can be calculated which comes from other elements. The numbers 
which represent this interaction field are preserved in the computer memory. 
After that an element is sought which is most wiling to flip. This element has 
to fulfil three conditions: 

1. Magnetization of this element must have opposite direction to the total 
field Htot exerted on this element (Htot = Hi + H ex ). 

2. Total field must have greater value than the switching field of the element 
(| H tot \> H s ). 

3. This element has maximal value of difference between the total field ex- 
erted on the element and its switching field. 

The selected element is fliped, and all numbers which represent the interaction 
fields are corrected by the change of the interaction field coming from the flipped 
element. 
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This procedure is repeated until no one element can be nipped. Then, the 
external field is changed. Subsequent changes of this field are of the values 
which make possible a flip of at least one element. The whole procedure is 
repeated. Finally the hysteresis loop of the whole array is obtained. 




Figure 1: The Pardavi-Horvath algorithm 



The algorithm is checked by the calculation of thehysteresis loop (see Fig.|5J 
and a comparison with the experimental one — the accordance is quanti- 
tative and good. 



3 Results of simulation 

3.1 The damage spreading 

The damage spreading technique is applied for an array of 100 x 100 elements. 
Two arrays are stored simultanously in the computer memory, in the same 
randomly selected initial state. The PH procedure is applied to lead these 
arrays to stable (stationary) states, where no flips occur. Then, one magnetic 
moment is flipped in one array, and wecheck if the flip is stable. Subsequently 
we apply a periodic magnetic field of amplitude H m . We consider the case of 
small frequency, where the whole arrays go to stable states each time before the 
field is varied. In this case a particular value of the frequency is not relevant. 
The damage is defined as the Hamming distance between two arrays, defined as 
the number of elements of the arrays with opposite directions of the magnetic 
moments. 
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Figure 2: The calculated hysteresis loop 



The result is that the Hamming distance increases only during some transient 
time. Then, the system reaches a limit cycle, with the length equal to some 
multiple of the period of the applied field. A typical result is shown in Fig. 
We note that if a system is critical, the size of damages should increase until 
the system boundaries are reached [HI El • In our case damages remain limited 
to a closed area at the lattice center. 

The size A of this maximal damage has a maximum for the amplitude of the 
applied field H m close to 1300 Oe, and it vanishes in most cases for H m > 1400 
Oe. It is obvious that any damage must disappearat strong fields, when the 
system is saturated. Besides, A varies strongly from a sample to asample. As 
a rule, the damage is localized as a shapeless and seemingly random pattern of 
the array elements, formed of several clusters, separated but close to each other. 

3.2 The fractal dimension 

The fractal dimension D of the cluster of the "up" (or "down" ) oriented elements 
is calculated for the 1024 x 1024 array. Obviously, the investigated magnetic state 
of the system should not be arbitrary, but it should be a consequence of the 
physical properties and processes. For a random initial state the fractal dimen- 
sion D is equal to 2.0, what means that the system has no fractal character. On 
the other hand, such a state is an artifact of the (pseudo)random numbers gen- 
erator, and therefore it cannot be treated as realistic. The only physical condition 
of such a state is its stability, but in our opinion this condition leaves too much 
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Figure 3: Time dependence of the damage. Maximal value of A in the stationary 
state is close to fifty. 



freedom. Therefore we have calculated the time evolution of the fractal dimen- 
sion D during the virtually performed hysteresis experiment. It is obvious, that 
D = 2.0 or D = 0.0 in the saturated states. We found that this bistable charac- 
ter of D is preserved also for small amplitude H m of the applied external field. 
Between these two values, the scaling is not proper. An example is shown in 
Fig. Ufa). Typically, a part of the curve for small size of the box shows the 
inclination close to zero, and another part — close to two. Similar plots have 
been presented also in Refs. [21 E] - As the applied field changes, the bent part 
of the curve is shifted left and finally disappears. In fact, plots like Fig. 0fa) 
appear when we try to assign the fractal character to a random configuration 
of spins. We demonstrate it in Fig. HJb), where the magnetic configuration of 
the array produces the curve shown in Fig. Ufa). Surely, this cluster does not 
exhibit any self-similarity. 

4 Conclusions 

The fractal dimension D obtained in Refs. [21 EI is 1.97, which is close to non- 
fractal value 2.0. On the other hand, in the Ising ferromagnet at the Curie point 
D = 187/96 w 1.95 [T3"] , Then, the value of D itself does not allow to state if 
the system is in the critical state. 

However, the criterion of criticality can be drawn from the obtained curve. 
Its bent character means that there is a characteristic length in the investi- 
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Figure 4: (a) The fractal dimension D as read from the slope log(n) vs log(d), 
where n is the number of boxes of size d, containing elements with "up" -oriented 
magnetic moments. Here D switches from two to zero, (b) A typical example of 
the array configuration for H = —400 Oe. The dots correspond to "up" -oriented 
array's elements. 



6 



gated structure. The value of this length is just at X coordinate of the bend 
point. We demonstrate it in Fig. [5] The two curves are obtained with the same 
box-counting technique for the conventional Ising ferromagnet. The curve for 
temperature close to the T c is a fairly straight line, except its upper part where 
a continuous deviation from linearity is observed. This deviation can be as- 
signed to a short waiting time (equivalent to N run ), which should be very large 
in the critical region, and to a finiteness of the system. The curve obtained for 
T = 0.44T C displays the same characteristic length, as in Fig. Ufa). This is an 
indication, that the system is not scale-free, and therefore not in a critical state. 
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Figure 5: The results of the box-counting technique to the Ising ferromagnet. 
Two curves refer to temperatures 1.01T C and 0.44T C . The straight line is a fit 
for T — 1.01T C . The result of the fitting gives D = 1.91, which is close to the 
theoretical value 1.95 



This conclusion is confirmed by the results obtained by means of the damage- 
spreading technique. In all cases investigated, we have not observed a damage 
which reaches the boundaries of the system. 

However, this conclusion is valid only for the system investigated experi- 
mentally in Ref. p^. We stress this point because there are serious differences 
between the this system and the model system discussed in Refs. |21 E| ■ In par- 
ticular, in the latter the lattice is triangular, what makes the system frustrated 
because of the antiferromagnetic-likc character of the magnetostatic interaction. 
Frustration enhances both the ground state energy and its degeneration. The 
system can wander in the phase space, and therefore damages can spread more 
easily. Then, the existence of SOC in an array with the triangular structure 
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remains to be checked. 
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